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ABSTRACT 


As  compared  to  the  commonly  used  constant  frequency  radar  waveforms,  the 
stepped-frequency  waveform  can  achieve  high  range  resolution  while  still  retaining  the 
advantages  of  lower  instantaneous  receiver  bandwidth  and  lower  analog-to-digital 
sampling  rate.  However,  the  relative  radial  motion  between  the  target  and  the  stepped- 
■frequency  radar  will  result  in  performance  degradations,  such  as  range  error,  loss  in 
signal-to-noise  ratio,  and  degraded  range  resolution.  The  solution  to  this  problem  is  to 
apply  velocity  compensation  to  the  received  signal,  which  can  eliminate  the  degradations 
due  to  Doppler  effects.  Three  velocity  compensation  schemes  for  the  detection  of  a 
moving  target  in  clutter  are  designed,  discussed,  and  compared  in  this  thesis.  Also,  a 
simulation  is  presented  to  verify  the  concepts,  and  simulation  results  are  compared  and 
discussed. 
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1.  INTRODUCTION 


The  stepped-frequency  waveform  is  an  interpulse  version  of  pulse  compression 
techniques  in  high  resolution  radars.  In  the  stepped-frequency  waveform,  the  carrier 
frequencies  for  a  burst  of  iV^  transmitted  pulses  are  monotonically  increased  from  pulse  to 
pulse  by  a  fixed  frequency  increment.  The  stepped-frequency  radar  can  achieve  high  range 
resolution  with  the  advantages  of  lower  instantaneous  receiver  bandwidth  and  lower 
analog-to-digital  (A/D)  sampling  rate,  both  of  which  are  the  major  limiting  factors  for 
achieving  high  range  resolution  in  the  conventional  constant  frequency  radar.  The  high 
range  resolution  of  the  stepped-frequency  radar  is  obtained  by  coherently  processing  a 
series  of  N  returns  from  the  stepped-frequency  pulses.  The  coherent  processing  is 
performed  by  taking  the  Discrete  Fourier  Transform  (DFT)  ofN  pulse  returns,  which  acts 
as  a  matched  filter  when  there  is  no  relative  radial  motion  between  the  radar  and  the 
target.  The  relative  radial  motion  during  the  coherent  processing  interval  will  result  in 
nonlinear  phase  shifts  in  the  received  signal.  As  a  result,  the  DFT  processing  is  not 
matched  to  the  received  signal,  leading  to  a  distorted  range  profile.  The  solution  to  this 
problem  is  to  apply  velocity  compensation  to  the  received  signal  based  on  the  target 
velocity.  This  thesis  mainly  focuses  on  the  study  of  velocity  compensation  in  the 
stepped-frequency  radar.  The  problem  of  the  detection  of  a  moving  target  in  clutter  will  be 
discussed,  which  requires  clutter  cancellation  to  eliminate  the  clutter  and  velocity 
compensation  to  restore  the  target  resolution. 


Three  velocity  compensation  schemes  are  designed  and  discussed  in  this  thesis.  In 
the  first  scheme,  the  velocity  compensation  is  performed  by  multiplying  by  a  compensation 
factor  in  the  time  domain,  while  in  the  second  scheme,  the  velocity  compensation  is 
performed  by  circular  convolution  with  the  DFT  form  of  the  compensation  factor  in  the 
fi'equency  domain.  Both  schemes  are  conceptually  equivalent,  but  have  different 
implementations.  Their  implementations  will  be  analyzed  and  compared.  In  the  first  and 
the  second  schemes,  the  target  velocity  is  assumed  to  be  known  in  advance.  The  third 
scheme  is  developed  based  on  the  first  scheme  with  the  addition  of  a  velocity 
compensation  loop  which  enables  the  target  with  unknown  velocity  to  be  processed; 
however,  the  target  velocity  is  assumed  within  a  defined  interval.  In  addition,  a  computer 
simulation  is  presented  to  evaluate  the  detection  of  a  moving  target  in  clutter  using  the 
stepped-fi'equency  waveform  by  means  of  the  three  velocity  compensation  schemes.  The 
simulation  results  are  discussed  and  compared. 

In  Chapter  H,  the  principle  of  the  stepped-ffequency  waveform  and  its 
corresponding  radar  system  are  introduced.  The  discussion  of  high  range  resolution 
processing  and  the  analysis  of  related  radar  parameters  are  also  included.  In  Chapter  III, 
three  velocity  compensation  schemes  for  the  detection  of  a  moving  target  in  clutter  are 
investigated.  The  simulation  results  of  the  three  compensation  schemes  are  also  presented 
in  this  chapter.  Conclusions  from  the  work  are  provided  in  Chapter  IV.  The  signal 
modeling  and  the  computer  programs  used  in  the  simulation  are  shown  in  Appendices  A 
and  B. 
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n.  PRINCIPLES  OF  THE  STEPPED  FREQUENCY  RADAR 


A.  STEPPED  FREQUENCY  WAVEFORM 

The  stepped-frequency  radar  transmits  a  burst  of  V  pulses,  whose  carrier  frequencies 

are  monotonically  increased  from  pulse  to  pulse  by  a  fixed  frequency  step  size,  as 
described  in  Figure  2.1.  Each  pulse  has  the  same  pulse  width,  x,  and  is  transmitted  at  a 
fixed  pulse  repetition  frequency  (PRF).  The  carrier  frequency  of  the  Arth  transmitted  pulse 
is  given  by 

where 

=  nominal  carrier  frequency 
A/"  =  frequency  step  size 
n  =1,  2,  ...,V 

The  frequency  spectrum  constituted  by  a  series  of  returns  from  a  target  due  to  the 
stepped-frequency  waveform  is  shown  in  Figure  2.2.  Thus,  the  total  effective  bandwidth  of 
the  stepped-frequency  waveform  is  given  by 

B^=NKf  (2.2) 
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Figure  2.1  The  Stepped-FrequencyWaveform. 


Figure  2.2  Total  Effective  Bandwidth  of  the  Stepped-Frequency 
Waveform. 
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The  S-dB  instantaneous  bandwith  of  the  stepped-frequency  waveform  is 
approximately  equal  to  the  inverse  of  the  pulse  width,  1/x,  and  is  much  less  than  the  total 
effective  bandwidth  of  the  waveform.  It  is  noted  that  the  range  resolution  of  the  stepped- 
frequency  waveform  radar  is  not  the  conventional  range  resolution,  ct/2.  The  synthetic 
range  resolution  achieved  by  coherently  processing  a  series  of  N  returns  from  a  target 
depends  upon  the  total  effective  bandwidth  and  is  given  as 


where 


c 

2NAf 


c  =  the  speed  of  light,  3x10*  m/sec 


(2.3) 


The  synthetic  range  resolution,  Ar,  given  by  Eq.2.3  for  a  specific  target  can  be 
improved  either  by  increasing  the  number  of  transmitted  pulses  in  a  burst,  N,  or  by 
increasing  the  frequency  step  size,  Af. 

B.  STEPPED  FREQUENCY  RADAR  SYSTEM 

The  implementation  of  the  stepped-frequency  radar  is  similar  to  that  of  a  coherent 
pulsed-Doppler  radar.  A  block  diagram  of  a  stepped-frequency  radar  system  is  shown  in 
Figure  2.3.  The  major  changes  in  the  stepped-frequency  radar  are  the  addition  of  a  fast 
coherent  stepped-frequency  synthesizer,  up-conversion  and  down-conversion  circuitry, 
and  the  requirement  for  a  wideband  fi’ont  end  [Ref  2:  p.236].  The  core  of  the  system  is  a 
coherent  stepped-frequency  synthesizer  with  the  output  frequency  of  The  frequency 
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of  the  synthesizer  is  stepped  from  pulse  to  pulse  by  a  fixed  frequency  step  size,  A/,  and  is 


written  as 


4.=("-W 

The  stepped-frequency  waveform  is  generated  by  first  mixing  the  frequency  of  a 
coherent  oscillator, with  the  output  frequency  of  the  stepped-frequency  synthesizer, 
f^.  Following  this  up-conversion,  the  signal  is  converted  to  the  final  transmitted  frequency 
by  mixing  with  the  frequency  of  a  stable  local  oscillator,  The  sum  frequency  signal  is 
then  pulse  modulated  and  amplified.  Thus,  the  carrier  frequenq?  of  each  transmitted  pulse 
is  composed  of  three  elements:  the  frequency  of  the  coherent  oscillator,  the  frequency 

of  the  stable  oscillator,  and  the  variable  stepped  frequency,  from  the  stepped- 
frequency  s)mthesizer.  The  carrier  frequency  of  the  Trth  transmitted  pulse  can  be  written  as 

+  (2-5) 

where 

nominal  carrier  frequency 

On  the  receiver  side,  the  return  signal  is  first  mixed  with  the  fixed  frequency  of  the 
stable  local  oscillator,/^,,  and  then  further  down-converted  by  mixing  with  a  sample  of 
the  stepped-frequency  S5mthesizer.  The  output  signal  after  the  second  mixer  is  the 


7 


intermediate-frequency  (IF)  signal  which  is  then  passed  through  an  IF  amplifier.  The  3-dB 

bandwidth  of  the  IF  amplifier  is  approximately  equal  to  the  inverse  of  the  pulse  width,  1/x, 

and  is  centered  on  the  frequency  of  the  coherent  oscillator,  The  stepped-frequency 

synthesizer  is  synchronized  with  the  transmitter  so  that  the  transmitted  carrier  frequency 

and  the  output  of  stepped-frequency  synthesizer  use  the  same  stepped  frequency, 

% 

within  a  pulse  repetition  interval  (PRI).  As  a  result,  multiple-time-around  target  echoes 
after  going  through  the  second  mixer  will  have  frequencies  varying  by  multiples  of  the 
frequency  step  size.  Therefore,  the  target  echo  whose  IF  signal  falls  outside  the  passband 
of  the  IF  amplifier  (i.e.,  target  echoes  from  different  transmitted  pulses  with  different 
will  be  filtered  out.  The  amplified  signal  is  then  split  into  two  channels,  the  in-phase 
channel,  I,  and  the  quadrature  channel,  Q,  which  are  mixed  with  the  frequency  of  the 
coherent  oscillator,  in  the  phase  detector.  The  in-phase  signal  is  mixed  with 
directly,  and  the  quadrature  signal  is  mixed  with  90°  shifted The  outputs  of  the  phase 
detector  are  sampled  and  quantized  in  the  A/D  converter.  It  should  be  noted  that  typically, 
the  sampling  interval  is  x  seconds  and  thus  each  sample  forms  a  range  cell.  The  samples 
from  the  A/D  converter  are  stored  in  memory  and  sorted  by  range  cells.  The  collected  N 
samples  from  a  range  cell  in  which  the  target  is  located  will  be  used  for  high  resolution 
processing  as  described  in  Figure  2.4. 
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Figure  2.4  N  Samples  from  a  Target  are  Collected  and  Stored  in 
Memory  for  High  Resolution  Processing. 

For  high  PRF  waveforms,  a  target  echo  may  arrive  in  a  different  PRI.  The  output  of 
the  stepped-frequency  synthesizer  will  have  switched  to  the  succeeding  stepped- 
frequencies  before  receiving  and  processing  the  previously  transmitted  pulses.  Two 
possible  solutions  to  this  problem  are  to  switch  the  synthesizer  to  the  previous  transmitted 
frequencies  or  to  have  several  synthesizers  to  sort  the  multiple-time-around  echoes  [Ref 


A  key  advantage  of  the  stepped-frequency  radar  over  both  the  ultra-wideband  radar 
and  the  typical  linear  frequency-modulated  (LFM)  pulse  compression  is  the  narrowband 
processing  which  eliminates  requirements  for  the  extremely  high  A/D  sampling  rate,  a 
major  constraint  for  conventional  pulse  compression  techniques.  The  disadvantage  of  the 
stepped  frequency  radar  is  that  it  takes  longer  time  to  transmit  and  receive  N  pulses  for 
high  resolution  processing.  However,  in  coherent  radars,  target  returns  from  #  transmitted 
pulses  are  collected  before  processing;  thus,  the  stepped-frequency  radar  does  not  pose 
longer  delay  than  the  normal  coherent  radar. 

C.  fflGH  RESOLUTION  PROCESSING 

The  sampled  output  of  the  phase  detector  for  the  /rth  pulse  can  be  written  as 

5(w)  =  a„  exp(/(J)„)  (2.6) 

where 

(J)„  =  the  phase  of  the  /rth  pulse 
a„  =  the  magnitude  of  the  /rth  pulse 

The  phase  component  of  the  sampled  output  for  the  /rth  pulse  is  given  by 

<|)„  =  27i/„  ^  (2.7) 

where  /,is  the  carrier  frequency  of  the  /rth  pulse  given  by  Eq.2. 1,  and  2R^/c  is  the  round 
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trip  time  of  the  target  return  due  to  the  wth  pulse.  R„  is  the  range  between  the  target  and 
the  radar  and  is  given  by 

R=R^-(n-l)vT  (2.8) 

where 

Rf^  =  initial  target  range 

T  =  pulse  repetition  interval 
V  =  relative  radial  velocity  between  target  and  radar 


The  expansion  of  the  phase  component  <()„  in  above  equation  contains  four  terms. 
The  first  term  of  Eq.2.9  is  the  constant  phase  resulting  from  the  initial  target  range  and  the 
nominal  carrier  frequency,  and  is  not  of  any  importance.  The  second  term  has  an  induced 
frequency  shift  by  the  amount  of  ILfRJTc,  resulting  from  the  frequency  step  size  during 
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the  round  trip  time  of  the  initial  target  range.  The  third  term  contains  a  Doppler  frequency 
shift  of  the  amount  Ic  due  to  the  relative  radial  velocity  between  the  radar  and  the 
target.  The  final  term  is  a  frequency  dispersion  which  is  determined  by  the  amount  of 
2v(w-l)A/c  due  to  the  relative  radial  velocity  and  the  output  frequency  of  the  frequency 
synthesizer.  The  third  and  fourth  terms  of  Eq.2.9  are  the  phase  shifts  caused  by  the 
relative  radial  velocity  between  the  radar  and  the  target,  which  will  result  in  target 
distortion  in  the  range  profiles. 

1.  Generation  of  HRR  profile 

After  target  samples  from  N  pulses  are  collected,  the  high  range  resolution  is 
obtained  by  coherently  processing  N  samples.  In  conventional  pulsed-Doppler  radars, 
coherent  processing  is  implemented  by  DFT  as  a  matched  filter.  Similarly,  the  high 
resolution  range  (HRR)  profile  in  the  stepped-frequency  radar  is  obtained  by  taking  the 
A-point  DFT  of  the  collected  A"  target  samples,  which  is  written  as 

S{k)  =  L  exp(/(j)„)exp  (-7^^) 

=  exp(,'2x/„2§i)exp(-y2^)  *=1,2,...,#  (2.10) 

Note  that  the  HRR  profile  maps  a  specific  range  cell  in  which  the  target  is  located. 
Thus,  the  absolute  range  of  the  target  is  not  shown  in  the  HRR  profile.  Only  the  relative 
positions  of  the  scatterers  on  the  target  are  determined.  The  DFT  bins  in  the  HRR  profile 
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can  be  interpreted  as  the  relative  synthetic  range  cells  which  resolve  the  relative  positions 
of  the  target  scatterers.  It  is  also  noted  that  the  same  procedure  for  obtaining  the  HRR 
profile  can  be  apphed  to  all  other  range  cells. 

To  reduce  the  sidelobe  effects  in  the  HRR  profile  caused  by  the  DFT  process,  the 
target  signal  is  multiplied  by  a  weighting  function  prior  to  computing  DFT.  In  this  thesis, 
Hamming  windows  are  employed.  The  weighting  function  can  suppress  the  range 
sidelobes  of  the  HRR  profile;  however,  it  also  causes  a  degradation  of  the  resolution  as 
indicated  by  a  scaling  factor  k  in  the  following  expression  for  the  range  resolution. 


where 


Ar  = 


kc 

2NAf 


it  =  1.33  for  Hamming  window 


(2.11) 


Note  that  the  DFT  process  functions  as  a  matched  filter  for  the  target  returns  due 
to  a  burst  of  N  transmitted  pulses  only  in  case  of  zero  relative  radial  velocity  between  the 
radar  and  the  target.  The  non-zero  relative  radial  velocity  will  result  in  a  mismatch 
between  the  DFT  and  the  Doppler-shifted  target  returns  [Ref  3:  p.4].  The  result  of  the 
mismatch  in  the  DFT  process  leads  to  target  spreading,  shifting  and  attenuation  in  the 
range  profile.  The  target  spreading  and  attenuation  are  considerably  more  important  than 
the  target  shifting  because  they  degrade  the  desired  range  resolution  and  decrease  the 
signal-to-noise  ratio.  The  effects  of  the  relative  radial  velocity  and  the  required  velocity 
compensation  will  be  discussed  in  the  next  section. 
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2.  Maximum  Unambiguous  Range, 

The  stepped-frequency  radar  achieves  high  range  resolution  by  generating  a  phase 
shift  by  means  of  linear  changes  in  the  carrier  frequencies.  The  induced  phase  shift  due  to 
the  frequency  change  from  pulse  to  pulse  is  given  by  the  second  term  of  Eq.2.9  as  follows 

(t.w  =2n^^(n-l)r  (2  12) 

The  induced  frequency  shift  is  obtained  from  above  equation  by  determining  the  rate  of  the 
phase  change  as 


(2.13) 

The  maximum  unambiguous  range,  is  obtained  by  substituting  the  maximum  induced 
frequency  shift  equal  to  PRF  as 


PRF=^2f=- 


(2.14) 


Therefore,  the  maximum  unambiguous  range,  in  the  HRR  profile  is  given  by 


(2.15) 
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The  fraction  of  the  maximum  unambiguous  range  window  occupied  by  an  original 
range  cell  AR  is  given  by 


M 

Ru 


ct/2 

dlAf 


=  tA/ 


(2.16) 


D.  EFFECTS  OF  THE  RELATIVE  RADIAL  VELOCITY 

The  third  and  fourth  terms  in  Eq.2.9  are  the  additional  phase  shifts  due  to  the 
relative  radial  velocity  between  the  radar  and  the  target,  resulting  in  a  mismatch  in  the 
DFT  process.  The  consequences  of  the  mismatch  in  the  DFT  are  three  negative  effects  on 
the  range  profile.  The  first  effect  is  target  spreading  in  range.  The  second  effect  is  the 
attenuation  of  target  amplitude  due  to  the  spreading.  The  amount  of  spreading  is 
determined  by  a  dimensionless  parameter  P  given  by  Eq.2.17  [Ref  2:  p.250],  which 
represents  the  number  of  the  synthetic  range  cells  (or  DFT  bins)  that  the  target  signal 
spreads  in  the  range  profile  during  the  coherent  processing  interval. 


P_vNT_  vN 
Ar  Ar(PRF) 


(2.17) 


The  spreading  of  the  target  signal  results  in  a  loss  of  synthetic  range  resolution,  for  the 
target  signal  has  broadened  from  a  synthetic  range  cell  with  resolution  Ar  up  to  several 
synthetic  range  cells. 
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The  third  effect  is  the  shifting  of  the  target  signal  in  the  range  profile.  The  target 
position  will  be  shifted  by  L  synthetic  range  cells,  where  L  is  given  by  Eq.2.18  [Ref.  3: 
p.45]. 


where 


L 


'an  a  a 

A  p_fcvNT_fc  vN 

Bejf  ^  5ej9-Ar(PRF) 


fc  =fo  +  ^Beff 


(2.18) 


If  the  number  of  the  shifted  synthetic  range  cells  L  is  greater  than  number  of  DFT  bins,  N, 
the  target  signal  in  the  range  profile  will  be  circularly  shifted,  resulting  in  a  "wraparound" 
[Ref  3:  p.41].  The  phenomena  of  the  target  spreading,  shifting  and  attenuation  in  the 
range  profile  caused  by  the  relative  radial  velocity  are  demonstrated  in  Figure  2.5. 

To  overcome  these  drawbacks  caused  by  the  relative  radial  velocity,  the  received 
signal  can  be  compensated  in  phase,  which  eliminates  the  effects  of  the  relative  radial 
velocity.  This  phase  compensation  depends  upon  the  target  velocity  and  is  normally 
termed  as  velocity  compensation  factor.  The  velocity  compensation  factor  will  consist  of 
the  third  and  fourth  terms  in  Eq.2.9  with  180  degrees  phase  rotation  as  given  by 

comp(n)  =  exp  (j^v(fo  +  (n-  l)A/)(«  -  l)r^  (219) 
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Magnitucle(dB) 


Phenomenons  of  target  spreading,  shifting,  attenuation 


Figure  2.5  Phenomena  of  Target  Spreading,  Shifting,  and 
Attenuation  for  a  Single  Target  Due  to  Differemt  Relative  Radial 
Velocities  in  the  Range  Profile.  The  Target  Range  is  Unchanged. 
SNR  =  13  dB,  PRF  =20  KHz,/„  =  IGHz,  N=  512,  A/=  1  MHz. 


With  the  application  of  the  velocity  compensation  factor  to  the  received  radar  signal,  the 
compensated  signal  becomes 

C(«)  =  5(«)exp  +  («  -  l)A/)(«  -  l)r)  (2.20) 

where  s(n)  is  given  by  Eq.2.6. 

For  a  moving  target  in  clutter,  with  DFT  processing,  the  moving  target  shifts 
from  the  clutter-filled  region  to  the  clutter-free  region  in  the  HRR  profile  as  shown  in 
Figure  2.6.  However,  the  straightforward  application  of  the  velocity  compensation  factor 
to  the  received  radar  signal  will  simply  compensate  the  target  signal,  but  the  clutter  will 
spread  in  the  range  profile  due  to  the  unnecessary  compensation  velocity.  Thus,  clutter 
cancellation  should  be  done  before  applying  the  velocity  compensation  factor  in  a  velocity 
compensation  scheme.  The  detailed  study  of  the  velocity  compensation  scheme  for  a 
moving  target  detection  in  clutter  using  the  stepped-frequency  waveform  will  be  discussed 
in  Chapter  HI. 

E.  SELECTION  OF  PULSE  WIDTH,  x  AND  FREQUENCY  STEP  SIZE,  A/ 

The  design  of  a  stepped-frequency  radar  system  requires  a  thorough  understanding 
of  the  effects  of  the  primary  waveform  parameters  on  the  system  performance.  The 
selection  of  the  primary  waveform  parameters,  such  as  the  pulse  width,  x  and  the 
frequency  step  size,  A^,  are  important  in  the  system  design.  The  selection  of  pulse  width,  x 
depends  on  the  maximum  target  extent,  E.  One  basic  requirement  for  good  target 
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Magnltude(dB) 


Figure  2.6  DFT  Processing  Makes  a  Moving  Target  Shift  to  the 
Clutter-Free  Region  in  the  Range  Profile.  The  Relative  Radial 
Velocity  =  500  m/sec,  PRF  =  20  KHz,  /„  =  lOGHz,  V=  512,  A/= 
1  MHz,  CNR  =  35  dB,  SNR  =11  dB. 
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delectability  is  that  the  pulse  width  must  encompass  the  maximum  target  extent.  This 


means  that  the  pulse  width,  x  must  be  greater  than  or  equal  to  the  round  trip  time  for 
length  E. 


T  >  ^  (2.21) 

or  E<^  (2.22) 

If  the  pulse  width  cannot  encompass  the  target  extent  (in  range),  different  scatterers  on  a 
target  will  not  be  located  within  the  same  original  range  cell,  resulting  in  a  loss  of  SNR  and 
an  incomplete  HRR  profile. 

The  requirement  in  choosing  the  frequency  step  size  Ef  is  that  the  maximum 
unambiguous  range  window  must  encompass  the  maximum  target  extent  E  to  avoid 
target  wraparound.  That  is, 


(2.23) 

or 

(2.24) 

When  dealing  with  moving  target  detection  in  clutter,  the  maximum  unambiguous  range 
window,  should  be  much  greater  than  the  original  range  cell,  NR  in  which  the  target  is 
contained.  That  is. 
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(2.25) 


or 


tA/‘«  1 


(2.26) 


Thus,  the  unused  portion  of  the  HRR  profile  is  a  clutter-fi-ee  region  which  can  be  used  to 
detect  the  shifted  target  signal.  In  this  situation,  the  original  range  cell,  AR  Avill  comprise 
only  a  small  fraction  of  the  maximum  unambiguous  range  window,  /?„. 

When  xA/  =  1,  the  maximum  unambiguous  range  window,  is  equal  to  an 
original  range  cell,  AR.  In  this  situation,  it  is  difficult  to  notice  the  target  wraparound 
caused  by  slight  changes  in  the  target  range.  Furthermore,  aliasing  may  occur  if  there  is 
interference  fi'om  adjacent  range  cells.  If  the  original  range  cell  encompasses  the  target  and 
the  clutter,  the  target  will  be  camouflaged  by  the  clutter  [Ref  1:  p.l68].  Thus,  the 
condition  of  xAf  =  1  is  only  effective  for  target  detection  in  a  non-clutter  environment. 

When  xAf  >  1,  the  maximum  unambiguous  range  window,  R^  is  smaller  than  the 
original  range  cell.  A/?.  As  a  result,  the  region  of  the  original  range  cell  not  contained 
within  the  maximum  unambiguous  range  window  will  fold  over  and  result  in  aliasing.  This 
situation  must  be  avoided  for  mapping  HRR  profiles  [Ref  6:  p.602-603]. 
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ffl.  VELOCITY  COMPENSATION  IN  THE  STEPPED  FREQUENCY  RADAR 


Velocity  compensation  plays  an  important  role  in  the  stepped-frequency  radar 
system,  for  the  resolution  in  the  range  profile  degrades  when  there  is  a  non-zero  relative 
radial  velocity  between  the  radar  and  the  target.  A  target  in  presence  of  clutter  can  be 
clearly  seen  only  when  sufficient  relative  radial  velocity  makes  the  target  shift  into  the 
clutter-fi'ee  region  in  the  range  profile.  With  DFT  processing,  the  moving  target  spreads  in 
range  and  attenuates  in  magnitude.  The  target  spreading  factor  P  given  by  Eq.2.17 
determines  the  number  of  the  synthetic  range  cells  occupied  by  the  moving  target  in  the 
range  profile.  The  effective  synthetic  resolution  for  a  moving  target  is  PAr  instead  of  the 
theoretical  resolution,  Ar.  In  addition  to  the  loss  of  range  resolution,  the  excessive  target 
shift  due  to  the  relative  radial  velocity  may  result  in  target  wraparound.  The  number  of 
synthetic  range  cells  that  a  moving  target  shifts  in  the  range  profile  is  equal  to  the  target 
shifting  factor,  L  given  by  Eq.2. 18. 

In  order  to  overcome  the  drawbacks  of  target  spreading  and  attenuation  due  to  the 
relative  radial  motion  between  the  radar  and  the  target,  the  velocity  compensation  factor 
given  by  Eq.2. 19  can  be  directly  applied  to  the  received  radar  signal.  However,  this  will 
simply  compensate  the  target  signal,  while  the  clutter  spreads  in  the  range  profiles  due  to 
the  unnecessary  compensation  velocity.  Thus,  in  a  velocity  compensation  scheme,  the 
clutter  cancellation  process  should  be  done  prior  to  applying  the  velocity  compensation 
factor. 
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Three  velocity  compensation  schemes  are  designed  and  discussed  in  this  chapter.  In 
the  first  and  the  second  schemes,  the  relative  radial  velocity  is  assumed  to  be  known  in 
advance.  The  third  scheme  deals  with  the  target  of  unknown  velocity,  but  the  velocity  is 
assumed  within  a  defined  interval.  Computer  simulations  for  the  three  compensation 
schemes  are  also  presented  in  this  chapter. 

A.  VELOCITY  COMPENSATION  SCHEME  1 

The  first  compensation  scheme  for  detection  of  moving  targets  in  clutter  is  shown  in 
Figure  3.1.  In  this  scheme,  the  received  return  signal  in  the  time  domain  is  first  windowed 
prior  to  taking  the  DFT.  This  DFT  process  generates  the  uncompensated  range  profile.  At 
this  stage,  the  clutter  occupies  N(xAf)  synthetic  range  cells,  which  is  canceled  by 
multiplying  by  a  rectangular  gating  function.  The  modified  range  profile  is  then 
transformed  back  to  the  time  domain  signal  by  an  inverse  DFT  (IDFT).  The  velocity 
compensation  factor  is  now  applied  to  the  clutter-canceled  target  signal  in  the  time 
domain.  After  velocity  compensation,  the  target  signal  is  windowed  and  the  DFT  is 
performed  once  again.  This  final  DFT  generates  a  clutter-free  and  velocity-compensated 
HRR  profile. 

Note  that  in  this  scheme,  two  weighting  functions  are  employed  to  reduce  the  range 
sidelobes  caused  by  the  DFT  process.  The  first  weighting  function  is  necessary,  and  the 
second  weighting  function  gives  a  minor  improvement  in  the  range  resolution. 
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Figure  3.1  Velocity  Compensation  Scheme  1. 


B.  VELOCITY  COMPENSATION  SCHEME  2 

In  this  section,  an  alternative  scheme  is  developed  for  velocity  compensation.  This 
new  scheme  is  based  on  frequency  convolution  [Ref.  4:  p.427]  which  is  shown  in  Figure 
3.2  and  written  as 


[xiin)  •  X2(w)]  =  i  IDFT[Xi(A:)  *  X2(k)]  (3.1) 


where  x^(n)  and  x/w)  are  finite  duration  time  domain  sequences  both  having  the  same 
length  V,  while  X^{k)  and  X^{k)  are  the  DFTs  of  xfri)  and  x^{n),  having  the  same  length  N. 
Note  that  the  convolution  in  the  fi-equency  domain  is  a  circular  convolution,  not  the 
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common  linear  convolution.  Thus,  the  result  of  frequency  convolution  will  also  be  a 


sequence  of  length  N. 


Figure  3.2  The  Property  of  Frequency  Convolution. 


In  Figure  3.1,  the  product  of  Xi(w)  and  x^{n)  is  replaced  by  the  circular  convolution 
of  X,(;t)  and  X.lk)  and  the  succeeding  IDFT,  where  x,(«)  represents  the  output  of  the 
IDFT  and  X2(n)  represent  the  velocity  compensation  factor.  Thus,  the  first  scheme  is 
expanded  as  shown  in  Figure  3.3.  The  purpose  of  this  substitution  is  to  cancel  out  the  two 
operations  of  IDFT  and  DFT  in  Figure  3.1,  resulting  in  a  modified  scheme  as  shown  in 
Figure  3.4.  The  cancellation  of  these  operations  is  expected  to  speed  up  the  computation. 
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Figure  3.3  The  Expansion  of  Velocity  Compensation  Scheme  1 
According  to  Frequency  Convolution. 
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Figure  3.4  Velocity  Compensation  Scheme  2. 


C.  ANALYSIS  AND  COMPARISON  OF  SCHEME  1  AND  SCHEME  2 

Both  velocity  compensation  schemes  are  functionally  equivalent.  In  structure,  the 
major  difference  between  the  two  schemes  is  that  in  scheme  1,  the  velocity  compensation 
is  done  by  multiplication  in  the  time  domain,  while  in  scheme  2,  the  velocity  compensation 
is  done  by  circular  convolution  in  the  frequency  domain.  As  a  basis  for  comparison,  we 
will  use  the  number  of  complex  multiplies  required  in  each  scheme  as  a  measure  of 
computational  complexity. 

In  scheme  1  (shown  in  Figure  3.1),  there  are  two  Hamming  windows,  three  DFT 
processes  (including  IDFT),  one  clutter  cancellation  filter,  and  one  complex  multiplication 
with  the  velocity  compensation  factor.  Assume  that  there  are  N  samples  for  the  receive 
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signal  and  AT is  a  power  of  2.  The  DFT  process  of  scheme  1  is  achieved  by  the  Fast  Fourier 
Transform  (FFT)  algorithm.  Thus,  there  are  N  complex  multiplies  in  each  Hamming 
window,  (N/2)\og^  N  complex  multiplies  in  each  DFT  process  [Ref  4;  p.506],  N  complex 
multiplies  in  the  clutter  cancellation  and  N  complex  multiplies  in  the  multiplication  with  the 
velocity  compensation  factor.  There  are  a  total  of  (3N/2)log^  N  +  AN  complex  multiplies 
involved  in  scheme  1  provided  that  an  FFT  algorithm  is  employed  in  the  DFT  process.  For 
example,  when  N  —  512,  there  will  be  a  total  of  8960  complex  multiplies  required  in 
scheme  1.  If  the  DFT  process  is  not  implemented  by  an  FFT  algorithm,  but  by  the 
standard  DFT  algorithm,  there  will  be  complex  multiplies  in  each  DFT  process  [Ref  4: 
p.426].  Thus,  there  will  be  a  total  o^3N'^  +  AN  complex  multiplies  involved  in  scheme  1. 
When  AT  =  512,  there  are  a  total  of  788,480  complex  multiplies  required  (approximately  88 
times  the  computational  cost  of  using  the  FFT  algorithm).  Therefore,  the  FFT  algorithm 
used  in  the  DFT  process  of  scheme  1  can  provide  much  higher  computational  speed  than 
the  standard  DFT  algorithm. 

In  scheme  2  (shown  in  Figure  3.4),  we  can  find  one  Hamming  window,  two  DFT 
processes,  one  clutter  cancellation  filter  and  one  circular  convolution  with  the  velocity 
compensation  factor  in  the  frequency  dommn.  A  direct  approach  for  the  circular 
convolution  process  requires  complex  multiplies  [Ref.  4:  p.426].  Assume  that  the  FFT 
algorithm  is  used  in  the  DFT  process  of  scheme  2.  Thus,  there  are  a  total  of  (N  log2  A^  + 
2A  +  AT  ^  complex  multiplies.  When  A^  =  512,  there  will  be  267,776  complex  multiplies 
required  in  scheme  2,  which  are  much  more  than  those  required  for  the  FFT  algorithm  in 
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scheme  1 .  If  the  standard  DFT  algorithm  is  employed  in  scheme  2,  there  will  be  a  total  of 
+  IN  complex  multiplies.  Thus,  for  N  =  512,  there  are  787,456  complex  multiplies 
required,  which  is  about  the  same  computational  complexity  as  using  the  standard  DFT 
algorithm  in  scheme  1.  As  expected,  for  both  velocity  compensation  schemes,  the  FFT 
algorithm  will  offer  less  computational  complexity  and  higher  processing  speed  than  the 
standard  DFT  algorithm. 

From  the  previous  analysis,  it  is  found  that  when  the  FFT  algorithm  is  employed  in 
the  DFT  process  of  both  schemes,  the  circular  convolution  will  be  the  major  factor  that 
results  in  the  much  greater  computational  complexity  in  scheme  2.  Therefore,  scheme  1 
mechanized  with  the  FFT  algorithm  is  recommended  for  velocity  compensation  in  the 
stepped-frequency  radar. 

D.  VELOCITY  COMPENSATION  SCHEME  3 

The  third  scheme  is  proposed  to  deal  with  the  target  of  unknown  velocity. 
However,  the  target  velocity  is  assumed  within  a  defined  interval.  Scheme  3,  as  shown  in 
Figure  3.5  is  implemented  by  using  the  model  of  scheme  1  or  scheme  2,  with  the  addition 
of  a  velocity  compensation  loop.  Scheme  1  is  employed  here  as  the  base  model  because  it 
has  much  higher  processing  speed  than  scheme  2.  In  scheme  3,  the  stepped-frequency 
radar  will  automatically  try  several  compensation  velocities  to  generate  compensated  HRR 
profiles.  These  compensated  HRR  profiles  will  be  stored  in  memory  and  compared  to  find 
the  one  with  the  best  resolution  and  signal-to-noise  ratio  (SNR).  It  should  be  noted  that  in 
scheme  3,  the  compensation  velocities  which  begin  at  a  defined  initial  value  are 
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Figure  3.5  Velocity  Compemsation  Scheme  3. 


31 


monotonically  increased  by  a  velocity  step  size,  Av  in  each  loop.  The  selection  of  velocity 
step  size,  Av  and  the  number  of  compensation  loops  will  depend  on  the  size  of  the  defined 
velocity  interval  and  the  waveform  parameters. 

E.  SIMULATION 

Scenario:  A  moving  target  approaches  a  surface-based  stepped-frequency  radar  at 
nose  aspect  as  shown  in  Figure  3.6.  The  radar  platform  and  the  clutter  are  assumed 
stationary.  The  target  with  four  dominant  scatterers  on  it,  is  8  m  long  at  a  distance  of 
11250  m  from  the  radar.  The  ranges  of  the  scatterers  on  the  target  are  11251  m  for 
scatterer  1  (RCS  =  0.6  m^ ),  11253  m  for  scatterer  2  {RCS  =  0.5  m^  ),  11255  m  for 
scatterer  3  (RCS  =  1  m^)  and  1 1257  m  for  scatterer  4  (RCS  =  0.7  m^).  The  relative  radial 
velocity  between  the  target  and  the  radar  is  initially  assumed  to  be  750m/sec,  and  will  be 
varied  in  the  simulation  for  the  purpose  of  testing  the  effects  of  the  relative  radial  velocity. 
The  clutter-to-noise  ratio  (CNR)  is  35  dB.  The  variance  of  the  thermal  noise  which  is 
modeled  as  white  Gaussian  noise  is  set  to  0.0002.  The  signal-to-noise  ratio  (SNR),  13  dB 
is  based  on  the  power  ratio  of  scatterer  3  (RCS  =  1  m^ )  to  the  thermal  noise. 

Waveform  parameters  used  in  the  simulation  are  as  follows. 


Number  of  pulses  per  burst  (N) :  512 

Nominal  carrier  frequency  (Q:  10  GHz 

Frequency  step  size  (A/):  1  MHz 

Pulse  repetition  frequency  (PRF):  200  KHz 

Pulse  width  (x):  0. 1  |j,sec 
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Figure  3.6  The  Pictorial  Description  of  the  Scenario  for  the 
Detection  of  a  Moving  Target  in  Clutter. 


Calculations  of  pertinent  quantities  are  as  follows. 
Total  effective  bandwidth: 


B^=mf=  500x1  MHz  =512  MHz 


Maximum  unambiguous  range  in  the  HRR  profile: 


T>  —  ^  - 

“  2A/' 


3  X  10*  m/ sec 


2  X  10^  Hz 
Original  range  cell  size: 


150  m 


_cx  _  (3  X  10*  m/sec)  x  (0.1  x  10~^sec)  _ 


m 


Synthetic  range  cell  size: 

c  Bu 


Ar  = 


Weff  2NAf  N 


=  0.293  m 
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The  fraction  ofi?„  occupied  by  AR: 

4^  =  -^^=tA/=  0.1 

Ru  dlAf 

The  number  of  synthetic  range  cells  within  the  fraction  of  AT?  in  the  HRR  profile: 
^  x#=  (tA/)  X  A^=  0.1  X  512  =  51.2  «  52 


Frequency  resolution  of  DFT; 

1  PRF  200  X  10^  Hz 


A^xPRI  N 


500 


=  390.625  Hz 


The  maximum  unambiguous  range,  R^  of  the  HRR  profile  and  its  corresponding 
maximum  induced  frequency,  PRF,  in  the  DFT  is  described  in  Figure  3.7.  Simulation 
programs  are  written  in  MATLAB  where  the  DFT  process  is  computed  by  using  a  radix-2 
FFT  algorithm  when  the  length  of  the  sequence  is  a  power  of  2  [Ref  5:  p.271]. 


Figure  3.7  The  Maximum  Unambiguous  Range,  R^  in  the  HRR 
Profile  is  Equivalent  to  the  Maximum  Induced  Frequency  Shift, 


PRF. 
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1.  Effects  of  The  Relative  Radial  Velocity 

To  demonstrate  effects  of  the  relative  radial  velocity  in  the  range  profiles  clearly,  the 
relative  radial  velocity  is  varied  in  each  uncompensated  range  profile.  The  relative  radial 
velocities  to  be  tested  in  this  simulation  are  0  m/sec,  350  m/sec,  650  m/sec,  950  m/sec  and 
1250  m/sec.  The  simulation  results  are  shown  in  Figure  3.8  through  Figure  3 . 12. 

In  Figure  3.8,  the  target  signal  is  in  clutter  and  cannot  be  detected  since  there  is  no 
relative  radial  velocity  between  the  target  and  the  radar.  In  Figure  3.9,  the  target  signal 
fi'om  a  moving  target  shifts  fi'om  the  clutter-filled  region  to  the  clutter-free  region,  but  the 
magnitude  of  target  signal  is  much  smaller  than  that  of  the  clutter.  Thus,  the  target  will 
possibly  be  treated  as  noise  and  neglected.  As  shown  in  Figure  3.10  through  3.12,  the 
scatterers  on  the  target  cannot  be  resolved,  and  the  phenomena  of  the  target  spreading, 
shifting,  and  signal  attenuation  in  the  range  profiles  become  more  prominent  as  the  relative 
radial  velocity  increases.  It  is  also  noted  that  in  Figure  3.9,  the  relative  radial  velocity  is 
350  m/sec,  but  the  scatterers  on  the  target  can  still  be  resolved  because  the  use  of  a  high 
PRF  (200  KHz)  reduces  target  spreading. 
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Figure  3.8  The  Uncompensated  HRR  Profile  for  a  Moving  Target 
with  Scatterers  in  Clutter.  The  Target  are  Masked  by  the  Clutter 
Since  There  is  No  Relative  Radial  Velocity  between  the  Radar  and 
the  Target. 


Figure  3.9  The  Uncompensated  HRR  Profile  for  a  Moving  Target 
with  Scatterers  in  Clutter.  The  Relative  Radial  Velocity  =  350 


m/sec. 


Magnitude(dB) 


50 


clutter 


Uncompensated  HRR  profile 


Figure  3.11  The  Uncompensated  HRR  Profile  for  a  Moving  Target 
with  Scatterers  in  Clutter.  The  Relative  Radial  Velocity  =  950 
m/sec. 
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Figure  3.12  The  Uncompensated  HRR  Profile  for  a  Moving  Target 
with  Scatterers  in  Clutter.  The  Relative  Radial  Velocity  =  1250 
m/sec. 

2.  Simulation  of  Scheme  1  And  Scheme  2 

In  this  subsection,  velocity  compensation  scheme  1  and  scheme  2  designed  for  the 
detection  of  a  moving  target  in  clutter  are  simulated.  The  radar  parameters  given  in  the 
scenario  are  used  in  the  simulation.  The  relative  radial  velocity  between  the  radar  and  the 
target  is  750  m/sec,  and  the  PRF  is  200  KHz.  Several  compensation  velocities,  such  as 
150  m/sec,  350  m/sec,  550  m/sec,  750  m/sec,  1000  m/sec,  and  1250  m/sec,  are  used  to 
evaluate  the  performance  of  both  compensation  schemes.  When  the  compensation  velocity 
is  equal  to  750  m/sec,  the  range  profile  is  correctly  compensated.  The  uncompensated 
range  profile  is  shown  in  Figure  3.13.  The  simulation  results  for  compensated  profiles  are 
shown  in  Figure  3.14  through  Figure  3.25. 
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In  Figure  3.14  through  Figure  3.25,  the  clutter  has  been  canceled  and  the  target 
resolution  becomes  finer  as  the  compensation  velocities  get  closer  to  the  relative  radial 
velocity,  750  m/sec.  In  the  correctly  compensated  HRR  profiles  shown  in  Figure  3.17  and 
Figure  3.22,  the  scatterers  on  the  target  are  the  most  clearly  resolved  and  have  the  best 
SNR,  for  the  relative  radial  velocity  is  fully  compensated,  and  hence,  the  DFT  process 
functions  as  a  matched  filter  to  the  received  signal.  Also,  the  SNR  (based  on  the  third 
target  scatterer  to  thermal  noise)  in  the  fully  compensated  HRR  profiles  is  almost  up  to  30 
dB  which  is  much  higher  than  the  previously  given  13  dB.  This  increased  SNR  is  due  to 
the  coherent  integration  of  ^pulses  in  the  DFT  process. 

Note  that  in  some  of  the  HRR  profiles  without  correct  velocity  compensation,  for 
example,  those  shown  in  Figures  3.16,  3.18 , 3.21,  and  3.23,  the  scatterers  on  the  target 
can  still  be  distinguished  because  the  PRF  used  in  the  simulation  has  a  high  value  of  200 
KHz  which  leads  to  lower  target  dispersion.  It  is  also  observed  fi’om  Figures  3.18  and  3.23 
that  the  target  signal  is  circularly  shifted  to  the  other  edge  in  the  range  profile  due  to  the 
excessive  compensation  velocity. 

In  general,  the  first  scheme  and  the  second  scheme  generate  very  similar 
compensated  HRR  profiles  for  the  same  compensated  velocity.  Nevertheless,  comparing 
the  simulation  results  of  both  schemes  carefully,  it  is  found  that  the  compensated  HRR 
profiles  generated  by  scheme  2  have  higher  range  sidelobes  and  noise  ripples.  The  lower 
range  sidelobes  in  the  compensated  HRR  profiles  generated  by  scheme  1  are  due  to  the 
existence  of  one  more  Hamming  window  in  the  scheme. 
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Figure  3.13  The  Uncompensated  HRR  Profile.  The  Relative  Radial 
Velocity  =  750  m/sec. 


Figure  3.14  The  Compensated  HRR  Profile  of  Scheme  1.  The 
Clutter  is  Canceled.  Compensation  Velocity  =150  m/sec. 
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Figure  3.15  The  Compensated  HRR  Profile  of  Scheme  1.  The 
Clutter  is  Canceled.  Compensation  Velocity  =  350m/sec. 


Figure  3.16  The  Compensated  HRR  Profile  of  Scheme  1.  The 
Clutter  is  Canceled.  Compensation  Velocity  =  550m/sec. 


Figure  3.17  The  Compensated  HRR  Profile  of  Scheme  1.  The 
Clutter  is  Canceled.  Compensation  Velocity  =  750  m/sec.  Note 
That  This  is  a  Correctly  Velocity-Compensated  HRR  Profile. 


Figure  3.18  The  Compensated  HRR  Profile  of  Scheme  1.  The 
Clutter  is  Canceled.  Compensation  Velocity  =  1000  m/sec.  The 
Target  Scatterers  are  Circularly  Shifted  to  Another  Edge  of  the 
Range  Profile  Due  to  the  Excessive  Compensation  Velocity. 
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Clutter  is  Canceled.  Compensation  Velocity  =  1250  m/sec.  The 
Target  Scatterers  are  Circularly  Shifted  to  Another  End  of  the 
Range  Profile  Due  to  the  Excessive  Compensation  Velocity. 


Figure  3.20  The  Compensated  HRR  Profile  of  Scheme  2.  The 
Clutter  is  Canceled.  Compensation  Velocity  =150  m/sec. 
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Figure  3.21  The  Compensated  HRR  Profile  of  Scheme  2.  The 
Clutter  is  Canceled.  Compensation  Velocity  =  350  m/sec. 


Figure  3.22  The  Compensated  HRR  Profile  of  Scheme  2.  The 


Clutter  is  Canceled.  Compensation  Velocity  =  550  m/sec. 
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Figure  3.23  The  Compensated  HRR  Profile  of  Scheme  2.  The 
Clutter  is  Canceled.  Compensation  Velocity  =  750  m/sec.  This  is  a 
Correctly  Velocity-Compensated  HRR  Profile. 


Relative  range(nieter) 


Figure  3.24  The  Compensated  HRR  Profile  of  Scheme  2.  The 
Clutter  is  Canceled.  Compensation  Velocity  =  1000  m/sec.  The 
Target  Scatterers  are  Circularly  Shifted  to  Another  End  of  the  HRR 
Proffle  Due  to  the  Excessive  Compensation  Velocity. 


Figure  3.25  The  Compensated  HRR  Profile  of  Scheme  2.  The 
Clutter  is  Canceled.  Compensation  Velocity  =  1250  m/sec.  The 
Target  Scatterers  are  Circularly  Shifted  to  Another  End  of  the  HRR 
Profile  Due  to  the  Excessive  Compensation  Velocity. 

3.  Simulation  of  Scheme  3 

In  the  simulation  of  scheme  3,  the  target  approaches  the  radar  platform  at  nose 
aspect  at  a  velocity  between  Mach  2  and  Mach  4  (approx.  660  m/sec  to  1320  m/sec).  The 
given  initial  compensation  velocity  is  660  m/sec,  and  the  velocity  step  size,  Av,  which  is 
added  in  each  compensation  loop  is  80  m/sec.  The  unknown  target  velocity,  which  is 
assumed  between  Mach  2  and  Mach  4,  is  generated  by  using  function  "rand"  in  MATLAB. 
The  simulation  results  are  shown  in  Figure  3.26  through  Figure  3.34. 

As  determined  from  the  simulation  results  of  compensation  scheme  3,  the 
compensated  HRR  profile  of  the  eighth  loop  (Figure  3.34)  has  the  best  target  resolution 
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Magnitude(dB) 


and  the  largest  SNR  Thus,  it  is  estimated  that  the  relative  radial  velocity  is  approximately 
equal  to  1220  m/sec. 


Figure  3.27  The  Compensated  HRR  Profile  of  The  First  Loop.  The 
Clutter  is  Canceled.  Compensation  Velocity  =  660  m/sec. 


Compensated  HRR  profile 


Figure  3.28  The  Compensated  HRR  Profile  of  the  Second  Loop. 
The  Clutter  is  Canceled.  Compensation  Velocity  =  740  m/sec. 


Magnitude(dB) 


Compensated  HRR  profile 


Figure  3  .29  The  Compensated  HRR  Profile  of  the  Third  Loop.  The 
Clutter  is  Canceled.  Compensation  Velocity  =  820  m/sec. 


The  Clutter  is  Canceled.  Compensation  Velocity  =  900  m/sec. 
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Figure  3.31  The  Compensated  HRR  Profile  of  the  Fifth  Loop.  The 
Clutter  is  Canceled.  Compensation  Velocity  =  980  m/sec. 


Figure  3.32  The  Compensated  HRR  Profile  of  the  Sixth  Loop.  The 
Clutter  is  Canceled.  Compensation  Velocity  =  1060  m/sec. 


Figure  3.33  The  Compensated  HRR  Profile  of  the  Seventh  Loop. 
The  Clutter  is  Canceled.  Compensation  Velocity  =  1 140  m/sec. 


Compensated  HRR  profile 


Figure  3.34  The  Compensated  HRR  Profile  of  the  Eighth  Loop. 
The  Clutter  is  Canceled.  Compensated  Velocity  =  1220  m/sec. 
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IV.  CONCLUSION 


In  the  stepped-frequency  radar,  the  high  range  resolution  is  achieved  by  performing 
the  DFT  on  the  received  target  signal.  Nevertheless,  this  is  optimum  only  when  there  is  no 
relative  radial  velocity  between  the  radar  and  the  target.  The  nonlinear  phase  shifts 
induced  by  the  relative  motion  cause  a  mismatch  in  the  DFT  process  and  result  in  a 
distorted  range  profile.  Direct  application  of  the  velocity  compensation  factor  to  the 
received  signal  allows  the  compensation  of  the  target  signal;  however,  the  clutter  signal 
will  disperse  in  the  range  profile  due  to  the  unnecessary  compensation  velocity.  Thus, 
clutter  cancellation  has  to  be  performed  before  velocity  compensation  in  a  velocity 
compensation  scheme.  In  this  thesis,  three  velocity  compensation  schemes  have  been 
presented  to  solve  this  problem. 

In  the  design  of  the  velocity  compensation  schemes,  weighting  functions  are  used  to 
suppress  range  sidelobes  in  the  compensated  range  profiles.  It  is  found  that  in  all  three 
schemes,  the  use  of  a  Hamming  window  prior  to  the  first  DFT  has  the  maximum  influence 
on  sidelobe  reduction.  In  addition,  the  second  compensation  scheme  was  originally 
expected  to  have  higher  computational  speed  than  the  first  scheme  due  to  fewer 
processing  stages.  However,  because  of  the  circular  convolution  of  the  compensation 
factor  and  the  target  signal,  the  complex  multiplies  used  in  the  second  scheme  are  much 
more  than  those  required  in  the  first  scheme,  resulting  in  significantly  longer  execution 
time.  The  first  and  the  second  schemes  are  theoretically  equivalent,  and  both  schemes 
generate  very  similar  compensated  HRR  profiles  for  the  same  compensation  velocity. 
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Nevertheless,  since  fewer  Hamming  windows  are  used  in  the  second  scheme,  it  is  found 
that  the  compensated  HRR  profiles  generated  by  the  second  scheme  have  higher  range 
sidelobes  and  noise  ripples.  These  phenomena  can  result  in  increased  false  alarms  or 
reduced  probability  of  target  detection.  From  the  viewpoint  of  fidelity  and  computational 
speed,  the  first  scheme  is  considered  better  than  the  second  scheme  for  implementation. 
The  third  scheme  is  based  on  the  first  scheme  with  the  addition  of  a  velocity  compensation 
loop.  In  this  case,  the  target  velocity  is  unknown,  but  is  assumed  within  a  defined  interval. 
In  addition,  it  is  observed  that  for  the  detection  of  moving  targets  using  the 
stepped-fi'equency  waveform,  the  use  of  a  high  PRF  is  appropriate  for  mapping  high  speed 
targets.  The  high  PRF  has  the  advantage  of  reducing  the  target  dispersion  even  if  accurate 
velocity  compensation  cannot  be  performed.  However,  a  high  PRF  mode  increases  the 
system  complexity. 

Waveform  design  for  the  stepped-frequency  radar  is  much  more  complex  than  that 
for  the  constant  fi'equency  radar.  The  selection  of  waveform  parameters  involves  several 
tradeoffs  that  are  not  present  in  standard  radars.  A  study  to  highlight  these  tradeoffs  and 
to  streamline  the  waveform  design  process  should  be  conducted  in  the  future. 
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APPENDIX  A.  SIGNAL  MODELING 


Target: 

In  the  scenario  of  the  simulation,  there  are  4  dominant  scatterers  on  the  target.  The 
magnitudes  of  N  return  pulses  from  a  specific  scatterer  on  the  target  are  assumed  constant. 
Thus,  the  sampled  scatterer  signal  of  the  phase  detector  for  the  «th  pulse  is  given  as 

Si  (n)  =  Ai  exp(/(j)„)  (A.  1) 

where 

Af  —  the  constant  magnitude  of  a  specific  scatterer 
for  the  wth  return  pulse 

(|)^  =  the  phase  component  of  a  scatterer  for  the  wth  pulse 
n  =  1,  2,  3, ...,  N 

i  -  index  i  is  used  to  distinguished  different 
scatterers  on  the  target,  i  =  1,  2,  3,  4 

The  phase  component  of  the  sampled  scatterer  signal  for  the  wth  pulse  is  given  by  Eq.2.7 
and  Eq.2.9  as  follows 


(j)„  =l%f„ 


2R„ 

c 


O  /y  /  .s  A  2(Ro  -(rt-  l)v7)^ 

—  27c(/o  +  (w  1)A^^  ^  J 
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where 


47c(/-o  +(rt-  l)AMRo  -(n-  l)v7) 
c 


(A.2) 


A/"  =  the  frequency  step  size  from  pulse  to  pulse 

V  =  the  relative  radial  velocity  between  the  target  and  radar 
(-  for  opening  velocity,  +  for  closing  velocity) 

T  =  pulse  repetition  frequency 

=  the  initial  range  when  target  is  detected 
c  =  the  speed  of  light,  3  x  1 0*  m/sec 

Thus,  a  sampled  scatterer  signal  is  modeled  as 


(«)  ^A,  exp 


(A3) 


Clutter: 

In  the  simulation,  it  is  assumed  that  the  stationary  clutter  distributes  in  an  original 
range  cell  in  which  a  moving  target  is  located.  Assume  that  there  are  15  clutter  sources, 
uniformly  distributing  in  an  original  range  cell  (AT?  =  15  m)  and  the  magnitudes  of  the 
return  pulses  from  all  15  clutter  source  are  constant.  Similarly,  the  sampled  clutter  signal 
of  each  clutter  source  for  the  wth  pulse  can  be  modeled  as 
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^  ( ATi{fo  +  (ri-\)^f)Rk\  /A /l^ 

Ck  («)  =  C  exp  ^ - c - ) 

where 

C  =  the  constant  magnitude  of  the  clutter  signal 

=  the  range  of  each  clutter  source 

k  =  index  k  is  used  to  distinguished  different 
clutter  sources,  A:  =  1,  2,  3, 15. 


Noise  Source  : 

The  thermal  noise  is  modeled  as  white  Gaussian  noise  with  variance  cP  The 
signal-to-noise  ratio  {SNR)  is  given  by 


SNR=^ 


(A.5) 


Thus,  the  constant  magnitude  of  a  scatterer  signal  can  be  wntten  as 

Ai  =  J2g^(SNR)  (A.6) 

The  clutter-to-noise  ratio  (CNR)  is  given  by 


CNR  = 


C^/2 


CL 

2a2 


(A.7) 
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The  constant  magnitude  of  the  clutter  signal  can  also  be  written  as 


C-  J^2a\CNR) 


(A.8) 


Overall  Received  Signal: 

The  overall  received  radar  signal  for  the  wth  return  pulse  from  an  original  range  cell 


is  modeled  as 


S(n)  =  E  Si  (n)  +  E  Cfe  (n)+noise(n) 


(A.9) 


where  noisein)  is  a  white  Gaussian  sequence. 
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APPENDIX  B.  SIMULATION  PROGRAMS 


function  [ru,  rstep]=schemel(N,fstep,prf,rl,r2,r3,r4,v) 

% 

%  This  program  is  a  MATLAB  function  which  computes  the  compensated  HRR  profile  of 
%  4-scatterer  moving  target  in  stationary  clutter  using  scheme  l.This  program  include 
%  clutter  cancellation  followed  by  velocity  compensation.  The  input  data  consists  of 
%  following: 

%  Number  of  pulses  (N) 

%  Frequency  step  size  (fstep,  Hz) 

%  Pulse  Repetition  Frequency  (prf,  Hz) 

%  Initial  Scatter  Ranges  (rl,r2,r3,r4,  meters) 

%  Target  Radial  Velocity  (v,  m/sec)(-opening)(+closing) 

%  The  following  parameters  can  be  changed  within  the  fimction 
%  Signal-to-Noise  Ratio  in  (dB)  (SNR) 

%  Clutter-to-Noise  Ratio  in  (dB)  (CNR) 

%  Nominal  Carrier  Frequency  (fo,  Hz) 

%  Velocity  Matched  for  Compensation  (co,  m/sec) 

%  The  outputs  are  Ru=c/2*fstep,  and  Rstep=Kc/2*n*fstep, 

%  which  are  the  maximum  unambiguous  range  and  processed  range  bin  size. 

%  The  FFT  size  is  equal  to  N. 

% 

% 


%  Radar  Parameters 

fo=10*10^9; 

n=0:N-l; 

c=3*10^8; 

f=fo+n*fstep; 

K=1.33; 

PRI=l/prf; 

cnr=35;  %dB 

snr=13;  %dB 

a=10^snr/10); 

A3=sqrt(2*a*0.0002); 

Al=sqrt(0.6)*A3; 

A2=sqrt(0.3)*A3; 

A4=sqrt(0.5)*A3; 

b=10^cnr/10); 

C=sqrt(2*b*0.0002); 
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%  Signal  Modelings,  Target,  Clutter  and  Noise 

sl=Al  *exp(lj*(4*pi/c)*f.  *(rH-n*v*PRI)); 

s2=A2*exp(lj*(4*pi/c)*f.*(r2+n*v*PRI)); 

s3=A3  *exp(  lj*(4*pi/c)*f.  *(r3+n*v*PRI)); 

s4=A4*exp(lj*(4*pi/c)*f.*(r4+n*v*PRI)); 

cl=C*exp(lj*(4*pi/c)*  1  *f); 

c2=C*exp(lj*(4*pi/c)*2*f); 

c3=C*exp(lj*(4*pi/c)*3*f); 

c4=C*exp(lj*(4*pi/c)*4*f); 

c5=C*exp(lj*(4*pi/c)*5*Q; 

c6=C*exp(lj*(4*pi/c)*6*Q; 

c7=C*exp(lj*(4*pi/c)*7*Q; 

c8=C*exp(lj*(4*pi/c)*8*Q; 

c9=C*exp(lj*(4*pi/c)*9*0; 

cl0=C*exp(lj*(4*pi/c)*10*f); 

cl  l=C*exp(lj*(4*pi/c)*  1 1  *f); 

cl2=C*exp(lj*(4*pi/c)*  12*f); 

cl3=C*exp(lj*(4*pi/c)*  1 3 

cl4=C*exp(lj*(4*pi/c)*14*Q, 

cl5=C*exp(lj*(4*pi/c)*15*Q; 

omegal=rand(size(n)); 

omega2=rand(size(n)); 

noise  =sqrt(0 .0002)*(omegal+lj*  omega2) ; 

s=s  1  +s2+s3  +s4+(c  1 +c2+c3 +c4+c5+c6+c7+c8+c9+c  1 0+c  1 1  +c  1 2+c  1 3 +c  1 4+c  1 5  )+noise; 
h=haniming(N)'; 

%  Uncompensated  HRR  profile 
s=s.*h; 

S=fft(s); 

ru=c/(2*fstep); 

rstep=(K*c)/(2*N*fstep) 

figure(l) 

plot(n,20*logl  0(abs(S))) 
grid 

titleCUncompensated  HRR  profile ') 
xlabelCRelative  range(meter)') 
ylabel('Magnitude(dB)') 
figure(2) 

t=linspace(l,ru,5 12); 

plot(t,20*logl0(abs(S))) 

grid 

title('Uncompensated  HRR  profile ') 
xlabel('Relative  range(meter)') 
ylabel('Magnitude(dB)') 
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%  Clutter  Cancellation 
hu=input('Enter  high  bin  number  for  clutter') 
co=input('Enter  compensation  velocity  for  target ') 
vc=exp(-lj*(4*pi*PRI/c)*(co)*f*n); 
ssl=[zeros(size(l:bu))  ones(size(l  :N-bu))].*S; 

%  Velocity  Compensation 
SSl=ifft(ssl).*vc; 

%  Second  Weighting 
SSl=SSl.*h; 

%  Compensated  HRR  profile 

FSS=fil(SSl); 

figure(3) 

plot(t,20*logl0(abs(FSS))) 

grid 

title('Compensated  HRR  profile') 
xlabelCRelative  range(meter)') 
ylabel('Magnitude(dB)') 
gtext('scatterers  1,2,3  and  4') 
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function  [ru,  rstep]=scheme2(N,fstep,prf,rl,r2,r3,r4,v) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

%  This  program  is  a  MATLAB  function  which  computes  the  compensated  HRR  profile  of 
%  4-scatterer  moving  target  in  stationary  clutter  using  scheme  2. This  program  include 
%  clutter  cancellation  followed  by  velocity  compensation.  The  input  data  consists  of 
%  following: 

%  Number  of  pulses  (N) 

%  Frequency  step  size  (fstep) 

%  Pulse  Repetition  Frequency  (prf) 

%  Initial  Scatter  Ranges  (rl,r2,r3,r4) 

%  Target  Radial  Velocity  (v)(-opening)(+closing) 

%  The  following  parameters  can  be  changed  within  the  function 
%  Signal-to-Noise  Ratio  in  (dB)  (SNR) 

%  Clutter-to-Noise  Ratio  in  (dB)  (CNR) 

%  Nominal  Carrier  Frequency  (fo) 

%  Velocity  Matched  for  Compensation  (co) 

%  The  outputs  are  Ru=c/2*fstep,  and  Rstep=Kc/2*n*fstep, 

%  which  are  the  maximum  unambiguous  range  and  processed  range  bin  size. 

%The  FFT  size  is  equal  to  N. 

% 

% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  Radar  Parameters 
fo=10*10^9; 
n=0:N-l; 
c=3*10^8; 

K=1.33; 

PRI=l/prf; 

cnr=35;%dB 

snr=13;%dB 

a=10^(snr/10); 

A3=sqrt(2*a*0.0002); 

Al=sqrt(0. 6)*  A3 ; 

A2=sqrt(0.3)*A3; 

A4=sqrt(0.5)*A3; 

b=10^(cnr/10); 

C=sqrt(2*b*0.0002); 

%  Signal  Modelings;  Targets,  Clutters  &  Noise 
f=fo+n*fstep; 

sl=Al*exp(lj*(4*pi/c)*£*(rl+n*v*PRI)); 

s2=A2*exp(lj*(4*pi/c)*£*(r2+n*v*PRI)); 


62 


s3=A3  *exp(lj*(4*pi/c)*f.  *(r3+n*v*PRI)); 
s4=A4*exp(lj  *(4*pi/c)*f.  *(r4+n*v*PRI)); 
cl=C*exp(lj*(4*pi/c)*  1  *f); 
c2=C*exp(lj*(4*pi/c)*2*f); 
c3=C*exp(lj*(4*pi/c)*3  *0; 
c4=C*exp(lj*(4*pi/c)*4*Q; 
c5=C*exp(lj*(4*pi/c)*5*Q; 
c6=C*exp(lj*(4*pi/c)*6*Q; 
c7=C*exp(lj*(4*pi/c)*7*Q; 
c8=C*exp(lj*(4*pi/c)*8*0; 
c9=C*exp(lj*(4*pi/c)*9*Q; 
clO=C*exp(lj*(4*pi/c)*  10*f); 
cl  l=C*exp(lj*(4*pi/c)*  1 1  *f); 
cl2=C*exp(lj*(4*pi/c)*  12*Q; 
cl3=C*exp(lj*(4*pi/c)*  13*f); 
cl4=C*exp(lj*(4*pi/c)*  14*f); 
cl5=C*exp(lj*(4*pi/c)*  15*^; 
omegal=rand(size(n)); 
omega2=Tand(size(n)); 
noise  =sqrt(0.0002)*(omegal+lj*omega2); 

S=sl+s2+s3+s4+(cl+c2+c3+c4+c5+c6+c7+c8+c9+cl0+cll+cl2+cl3+cl4+cl5)+noise; 

h=haniming(N)'; 

%  Uncompensated  HRR  profile 
s=s.*h; 

S=ffi(s); 

ru=c/(2*fstep); 

rstep=(K*c)/(2*N*fstep) 

figure(l) 

t=linspace(  1  ,ru,N); 
plot(t,20*logl  0(abs(S))) 
grid 

titleCUncompensated  HRR  Profile') 
xlabelCRelative  range(meter)') 
ylabel('Magnitude(dB)') 
figure(2) 

plot(t,20*logl  0(abs(S))) 
grid 

titleCUncompensated  HRR  Profile') 
xlabelCRelative  range(meter)') 
ylabelCMagnitude(dB)') 

%  Clutter  Cancellation 

bu=input(Enter  high  bin  number  for  clutter') 

co=input('Enter  compensation  velocity  for  target ') 
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vc=exp(- lj*(4*pi*PRI/c)*co*f.  *n); 
ss  l=[zeros(size(  1  :bu))  ones(size(  1  ;N-bu))] .  *  S; 
%  DFT  of  Velocity  Compensation  Factor 
vc=(  1/N)*ffl(vc); 

%  Velocity  Compensation 
SSI  =circonv(ssl ,  vc); 

%  Compensated  HRR  profile 
figure(3) 

plot(t,20*log  1 0(abs(S  SI))) 
grid 

title('Compensated  HRR  profile') 
xlabelCRelative  range(meter)') 
ylabel('Magnitude(dB)') 
gtext('scatterers  1,2,3  and  4') 
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function  [value]=circonv(xl,x2) 

% 

%  This  MATLAB  function  is  to  calculate  cicular  convolution  via  direct  approach.  Both  xl 
%  and  x2  are  sequences  of  length  N.  The  resulting  sequence  of  cicular 
%  convolution  has  also  the  same  length  N. 

% 

% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


k=2;length(x2); 

A(l,l)=x2(l); 

A(l,k)=x2(length(x2)-(k-2)); 
vjdue(l)=sum(xl.*A(l,:)); 
for  i=2:length(x2); 
for  k=2:length(x2); 

A(i,  l)=A(i-l,length(x2)); 
A(i,k)=A(i-l,k-l); 
end 

value(i)=sum(xl .  *  A(i, ;)); 
end 
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function  [ru,  rstep]=scheme3(N,fstep,prf,rl,r2,r3,r4) 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
% 

%  This  program  is  a  MATLAB  fiinction  which  computes  the  compensated  HRR  profile  of 
%  4-scatterer  moving  target  in  stationary  clutter  using  scheme  3. This  program  include 
%  clutter  cancellation  followed  by  velocity  compensation.  The  input  data  consists  of 
%  following: 

%  Number  of  pulses  (N) 

%  Frequency  step  size  (fstep) 

%  Pulse  Repetition  Frequency  (prf) 

%  Initial  Scatter  Ranges  (rl,r2,r3,r4) 

%  The  following  parameters  can  be  changed  within  the  function 
%  Signal-to-Noise  Ratio  in  (dB)  (SNR) 

%  Clutter-to-Noise  Ratio  in  (dB)  (CNR) 

%  Nominal  Carrier  Frequency  (fo) 

%  Initial  Compensation  Velocity  (v)(-opening)(+closing) 

%  Velocity  Matched  for  Compensation  (co) 

%  The  outputs  are  Ru=c/2*fstep,  and  Rstep=Kc/2*n*fstep, 

%  which  are  the  maximum  unambiguous  range  and  processed  range  bin  size. 

%  The  FFT  size  is  equal  to  N. 

% 

% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  Radar  Parameters 

v=-(660+660*rand(l)); 

fo=10*10^9; 

n=0:N-l; 

c=3*10^8; 

K=1.33; 

PRI=l/prf; 

cnr=35; 

snr=13; 

a=10^(snr/10); 

A3=sqrt(2*a*0.0002); 

Al=sqrt(0.6)*A3; 

A2=sqrt(0.3)*A3 ; 

A4=sqrt(0.5)*A3; 

b=10^(cnr/10); 

C=sqrt(2*b*0.0002); 

%  Signal  Modelings:  Targets,  Clutter  &  Noise 
f=fo+n*fstep; 

sl=Al*exp(lj*(4*pi/c)*f  *(rl-n*v*PRI)); 
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s2=A2*exp(lj  *(4*pi/c)*f.  *(r2-n*v*PRI)); 
s3=A3  *exp(  Ij  *(4*pi/c)*f.  *(r3  -n*v*PRI)); 
s4=A2*exp(lj  *(4*pi/c)*f.  *(r4-n*v*PRI)); 
cl=C*exp(lj*(4*pi/c)*  1  *0; 
c2=C*exp(lj*(4*pi/c)*2*f); 
c3=C*exp(lj*(4*pi/c)*3  *0; 
c4=C*exp(lj*(4*pi/c)*4*Q; 
c5=C*exp(lj*(4*pyc)*5*f); 
c6=C*exp(lj*(4*pi/c)*6*Q; 
c7=C*exp(lj*(4*pi/c)*7*f); 
c8=C*exp(lj*(4*pi/c)*8*0; 
c9=C*exp(lj*(4*pi/c)*9*0; 
clO=C*exp(lj*(4*pi/c)*  10*f); 
cl  l=C*exp(lj*(4*pi/c)*  11*0; 
cl2=C*exp(lj*(4*pi/c)*  12*f); 
cl3=C*exp(lj*(4*pi/c)*  13*f); 
cl4=C*exp(lj*(4*pi/c)*  14*0; 
cl5=C*exp(lj*(4*pi/c)*  15*0; 

omegal=rand(size(n)); 
omega2=Tand(size(n)); 
noise  =sqrt(0.0002)*(omegal+lj*omega2); 

S=sl+s2+s3+s4+(cl+c2+c3+c4+c5+c6+c7+c8+c9+cl0+cin-cl2+cl3+cl4+cl5)+noise 

h=hamining(N)'; 

%  Uncompensated  HRR  profile 
s=s.*h; 

S=fft(s); 

mag=abs(S); 

ru=c/(2*fetep); 

rstep=(K*c)/(2*N*fstep) 

figure(l) 

plot(n,20  *log  1 0(mag)) 
grid 

titleCUncompensated  HRR  Profile') 
xlabelCRelative  range(meter)') 
ylabel('Magnitude(dB)') 
figure(2) 

t=linspace(l,ru,5 12); 

plot(t,20*logl0(mag)) 

grid 

titleCUncompensated  HRR  Profile') 
xlabelCRelative  range(meter)') 
ylabel('Magnitude(dB)') 

%  Clutter  Cancellation 


bu^nputCEnter  high  bin  number  for  clutter') 
co=  -660; 

M=l:8; 

cco=co-(M-l)*80; 

forP=l:8 

vc=exp(lj*(4*pi*PRI/c)*cco(P)*£*n); 
ssl=[zeros(size(l :bu))  ones(size(l ;N-bu))]. *S; 
%  Velocity  Compensation 
SSl=ifB;(ssl).*vc; 

%  Second  Weighting 
SSl=SSl.*h; 

%  Compensated  HRR  profile 

FSS=ffl(SSl); 

figure(P+2) 

plot(t,20*logl  0(abs(FS  S))) 
grid 

title('Compensated  HRR  profile') 
xlabel('Relative  range(meter)') 
ylabel('Magnitude(dB)') 
end 
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